1 



Self-dual Ginzburg-Landau vortices in a disk 

G.S.Lozano a - fc , M.VMamas c ' fc and E.F.Moreno c < & 
a Depto. de Fisica, FCEyN, Universidad de Buenos Aires, Pab I, Ciudad Univeristaria, Buenos Aires, Argentina. 
■ b Consejo Nacional de Investigaciones Cientificas y Tecnicas, Argentina. 

c Depto. de Fisica. Universidad Nacional de La Plata. CC 67, 1900 La Plata, Argentina. 

(May 2000) 



(N 



a ■ 

We study the properties of the Ginzburg-Laundau model in the self-dual point for a two- 
dimensional finite system . By a numerical calculation we analyze the solutions of the Euler-Lagrange 
equations for a cylindrically symmetric ansatz. We also study the self-dual equations for this case. 
We find that the minimal energy configurations are not given by the Bogomol'nyi equations but by 
solutions to the Euler Lagrange ones. With a simple approximation scheme we reproduce the result 
of the numerical calculation. 
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The study of vortex solutions in Ginzburg Landau (GL) theories has been the subject of continuous interest in 
j3 ' different areas of condensed matter and high energy physics. 

Static solutions in 2-dimensional infinite samples are characterized by the dimensionless GL parameter k which is 
' . defined as the ratio of the magnetic penetration length A and the coherence length £. It has been known for several 
years that the GL model in the infinite plane possesses very special properties at n 2 — 1/2. For instance the second 
order static equations of motion are equivalent to a set of much simpler first order equations referred as self-dual or 
q ' Bogomol'nyi equations (BE) [Q. 

Although the existence of self-dual equations was first pointed out within the study of superconductors by Harden 
and Arp |2j most of the research on the subject has been done in connection with the role this type of equations play 
in High Energy Physics. Indeed, self-dual equations were introduced in this context independently by Bogomol'nyi 
and de Vega and Schaposnik jjj in their study of vortex solutions of the relativistic version of the GL model (Abelian 
^\ " Higgs model). Since then, the properties of the solutions, the connection with Supersymmetry, Topological Field 
t-H , Theories and Duality have been established not only for the Abelian Higgs model but also for related theories in 
*P ' different number of space time dimensions and non- Abelian gauge groups Q| . 

Very recently, Akkermans and Mallick have addressed the study of the GL model at the self-dual point for a 2- 
dimensional disk of finite radius R j|. Their study is relevant for the case of vortices in mesoscopic systems, where the 
size of the sample is of the order of A and As evidenced by recent experiments @] , the superconducting behavior 
of mesoscopic disks is radically different from that of the same material in the macroscopic regime || . An interesting 
question that arises is then to determine in which way size effects manifest at the self-dual point and reciprocally to 
analyze if the special properties that the model show at the self-dual point allow for a simpler interpretation of the 
'■^ ' experimental results 

In this work we shall re-analyze the properties of the GL model for finite systems in the self-dual point. We will 
perform a numerical study of the equation of motion and explore the role played by the self-dual equations for this 
case. As a result of our analysis it will be shown that some of the approximations made in Q are not correct. We 
shall in turn present a simple approximation scheme which correctly reproduces the numerical calculation. 
The GL expression for the energy of a 2-dimcnsional sample Q can be written as 
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E = J dx { + jWf+^WD) (!) 

where i*y = diAj — djAi, and Di<p = di<p — iqAi<p with i = 1,2. Here Ai denotes the electromagnetic vector potential, 
is a complex scalar field (order parameter) and q, the charge. Writing the potential as 

V(\<P\)=P/2{\^- V 2 ) 2 , (2) 

the penetration length and the coherence length are A 2 = \/(Airq 2 ri 2 ) and C 2 = 1 / (27/ 2 /3) while the GL parameter is 
K 2 =f3/(2nq 2 ). 

For arbitrary k the minimal energy configurations satisfy the second order Euler Lagrange equations, 
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DM = (3) 



^diFij = -jj = -q/2i{4>*D 4> - <j>D^*) (4) 



Using the identity 



\ D,o | ><,,/)/:■ - = \\D^\ 2 ± ^jdiJj ± |-B|.r| (■-.) 
the energy at the self-dual point can be re-written as 



where $ = f n d 2 x B = <f> dn A ■ dx is the total magnetic flux through the sample. Assuming that the current is zero at 
the boundary, a lower bound for the energy is obtained, 

E>\&*\. (7) 

Energy configurations saturating the bound must satisfy the self-dual or Bogomol'nyi equations, 

F ij ±2Trqe ij (\<t>\ 2 -r) 2 ) = , (8) 

Di<t> =F ieijDj(j) = 0. (9) 

In the plane these equations are totally equivalent to the Euler Lagrange equations. To obtain finite energy 
conditions, one has to demand 



which in turns implies that 



IimA$ = lim |$| 2 = r? 2 (fO) 

p— >OG p — >OC 



lim Ji = 0; (f f) 

p— >oo 



On the other hand, writing <p — |</>|e 2X , the current takes the form: 

J l =q\<P\ 2 (d l X-qA i ). (12) 

As is a single valued field, the phase x{pi 6) must satisfy 

X(p,27r)- X (p,0) = 27rn, (13) 

with n an integer. This condition together with (^) implies that the total flux has to be an integer multiple of the 
quantum of the flux $o = 2k /q 

$ = / Aidx 1 = - [ d lX dx i = n§ , (14) 
Jan q Jon 

For each integer n, equations (||), admit a family of solutions depending on 2\n\ parameters that can be identified 
as the 2 dimensional coordinates of |n| non interacting vortices with flux quantum $ (the upper or lower sign in the 
equations has to be chosen according to the sign of n)( @). 

Let us now concentrate to the case of finite geometries with a boundary. When the current is zero at the boundary, 
the Bogomol'nyi inequality still provides a bound for the energy and the flux is quantized in units of <!>o- However, 
for a given flux, there is a minimal area of the sample for which the equations admit solutions. Indeed, integrating 
the first equations we get the inequality: 
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$= / B= / 27rq(?7 2 - |0| 2 ) < / 2irr] 2 = ?? 2 | area(ft). (15) 
Jn Jn Jn 2 

In the infinite plane the requirement of zero current at the boundary is the natural boundary condition since it is 
the only way of obtaining finite energy solutions. However there is no compelling reason to do so for a finite region. 
The appropriate boundary condition is Q , 

D ± (j) = 0. (16) 

Notice that this boundary condition only implies vanishing of the normal component of the current at the boundary. 
The tangential component is left in principle undetermined. Nevertheless, it is possible to show that for configurations 
satisfying the self-dual equations, it holds the relation, 

J i = ±le ij d j \<f> 2 \ (17) 

That is, for solutions of the BE equations both components of the current must vanish at the boundary. As discussed 
above, this implies that the total flux is quantized. 

Suppose now that instead of imposing a boundary condition over Jn, we fix the total flux $. It is clear, that if 
$/$o is not an integer, a tangential component of the current will be established at the boundary and minimal energy 
configurations will not be given by solutions to the BE but by solutions to the Euler Lagrange equations. More 
interestingly, when <f> /$ — m i with m integer, although the BE do admit solutions, they are not the minimal energy 
configuration. In fact, it is energetically favorable to create n < m vortices and to develop a tangential current at the 
boundary. 

Let us illustrate this in a cylindrically symmetric ansatz: 

<t>(x) = f(f)e me 
A e {x)=A{p) 

A p (x) = . (18) 
Defining dimensionless variables x(r) = n — qA(r), z(r) = f(r)/r), r = (p/A), the self-dual equations become 

x'± r -{z 2 -l)=0 (19) 

z' T — = (20) 
r 

and the boundary condition ( fuj ) translates into a Neumann condition for the order parameter: 

z'(R*) = (21) 

where R* = R/X. Thus, we see that unless x(R*) = (which implies that $ = § A — 2im/q) the BE ( po| ) cannot be 
satisfied at r — R*. 

Minimal energy solutions are then obtained by solving the second order Euler-Lagrange equations. In our ansatz 
they read, 



d 2 x 1 dx 2 () 

dr 2 r dr 



d z 1 dz X A, ■■ , 

d^+rTr- — + KZ{1 - Z) = °> (23) 



while the energy can be expressed as, 



E = n V 2 I" rdr[{-) 2 + z' 2 + ^ + n 2 (z 2 - if]. (24) 



r r 2 



Regularity of cylindrical coordinates impose conditions at r = 0, and together with the Neumann condition ( ^l| ) and 
the definition of x(R*) in terms of n and $ we have the following set of boundary conditions: 
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x(0)=n x{R*)=n - — 

z(0)=0 z'(R*)=0. (25) 

We have analyzed the existence of solutions of the system (|2^,^3|,^5|) by numerical integration. We employed a 
relaxation method for boundary value problems In such method, the differential equations are discretized in a 
convenient mesh and converted to a set of coupled algebraic equations. The system is then solved using Newton's 
iterative method, starting from an initial guess and improving it iteratively. Given $ and R* a solution for each 
integer n is found corresponding to a local minimum of the energy. Then, the n giving the lowest energy is selected. 

In figure 1 we show the solutions corresponding to $ = 6$o> R* — 10- For these values, the minimal energy 
configuration corresponds to n = 1. The energy of this solution is E = 3.5l7r?y 2 . As the external flux is an integer 
number of flux quanta, the Bogomoln'nyi equations also have solutions for this case with energy E = 6nr] 2 . Thus, the 
system lowers its energy by allowing a tangential component of the current Jn(R*) ^ 0. 




FIG. 1. Vortex solution and angular component the current j(r) for R = 10, n = 1 and $ = 6$o- 



Notice nevertheless that there is a point Rq such that J(Rq) = 0. Following reference |Q, we will consider the disk 
as 51 = Oi U where Oi, is the inner disk < r < Rq and f^, the outer ring Rq < r < R*. We will express the 
energy as: 

E(Q) = E(Clx) + E(Cl 2 ) (26) 

As the current vanishes at Rq, the authors in Q assumed that the minimal energy configurations satisfy the 
Bogomol'nyi equations in f2i; however it can be shown that this is not the case. Consider the function L[r] = 
z'(r) =p x(r) z(r)/r, evaluated with the solutions of the Euler-Lagrange equations. If the Bogomol'nyi equations are 
satisfied on S7i , then L[r] = in this region. But, as the current is not zero on the external boundary r = R* , the 
Bogomol'nyi equations are not satisfied on f^, and thus L[r] ^ on ^2. Clearly, regularity of the solutions of ordinary 
differential equations prevents the existence of a function that vanishes in the whole region < r < Rq but is different 
from zero for r > Rq. Although the solutions of the self-dual equation minimize the energy on the internal region f2i, 
any regular extension of the solution to the whole disk will not minimize the total energy. 

Having said this we should notice that even though the self-dual solutions are not exact solutions in the inner disk 
fii, they are in fact a very good approximation. A numerical analysis of both solutions shows that for a wide range of 
parameters, they differ only in around one part in a thousand and the same is true for the energy. Then we will take, 

£(fi x ) w mf\n\ (27) 

Let us now analyze the contribution of the O2 region to the energy. 

E(Sh) = W £ rdr[{^f + z' 2 + ^ + \{z 2 - I) 2 ] (28) 

We first review the main steps in Qj. There it was assumed that as the fields are concentrated in a region of width 
of order one from the border, this expression could be approximated as 
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E(Q 2 ) « W r[C-f + z' 2 + ^- + ^(z 2 1) 



r=R* 



Then, the condition ^M- = would give 



x 2 (R) 
R 2 



\-z l 



(29) 



(30) 



As a next step Akkcrmans and Mallick neglected the magnetic energy contribution (first term in (|29|)) arriving to the 
following expression for the energy, 



E(n 2 ) ~ 7TTJ R 



2 ^ f x 2 (R*) l* 4 (i?*) 



R* 2 2 R* 4 ' 

Finally, the authors neglected, in the large R limit, the quartic term in ( |3l"| ) ending with the expression 



E(Q 2 ) « mj 



■ x 2 (R*) 
R* 



(31) 



(32) 



The reasoning above suffers from two main drawbacks. Although the field z is practically constant, the field x is 
not, making the approximation of the integral not valid. In fact, our numerical simulation shows that equation ( |30| ) 
is not fulfilled. Second, as shown in fig 2 the magnetic energy is of the same order than the third term in (Bq). 




FIG. 2. Energy density contributions near the edge for various terms (R* = 10, n = 1 and $ = 6$o)- The filled line 
corresponds to the term B 2 , the dashed-dotted line corresponds to the |0| 2 |Vx — A\ 2 term , and the dashed line to the term 
| \7<3> | 2 . Note that the B 2 term provides the largest contribution to the energy and thus cannot be neglected. 

The correct field distribution in £l 2 can instead be obtained by a simple approximation using 



R* 



(33) 



as an expansion parameter. Indeed, defining 



X = x I V 



(34) 



the equations of motion become 



d 2 x 1 dx 



dr 2 



dr 



z 2 x = 



(35) 



d z 
dr 2 



1 dz 
r dr 



+ 2^( 1 ^ 2 ) 



(36) 
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with the boundary conditions 



x{R* ) = x(R*)=R* 



(37) 



z{R* ) = 1 z'(R*) = 

In order to solve these equations, we make an expansion of the form 



x = x + v 2 x 2 + 0(v A ) 
z=l + v 2 z 2 + 0(v 4 ) 



(38) 



(39) 



The boundary conditions will be satisfied by the solution x while homogeneous conditions are valid for x 2 and z 2 . 
To lowest order we obtain 



d 2 x _ 1 cEo _ - _ Q 
dr 2 r dr 



with solution 



h{r)-^-Mr) 



Ki(R* ) 



(40) 



(41) 



^^Ik^R*)} and h{r), K^r) are Besscl functions. 



Here Cq 1 = [h{R*) ^ 

With these results we can solve the order v 2 equations: 



d 2 x 2 1 dx 2 
dr 2 r dr 



x 2 = 2z 2 x , 



(42) 



d 2 z 2 1 dz^ ^ 2 



2 z 2 = ^ (43) 



dr 2 r dr 

with homogeneous boundary conditions for both functions. Using the Green functions for each equation we obtain 
the solution 

Z2{r)=y 2 {r) — — dr + yi {r) — — rfr (44) 



where 



w(r)=^o(r)-^^/o(r) 
7 o(-ft ) 

y2(r)=Ko(r)-^^Io(r) (45) 



/(r) 



and W yi , V2 is the Wronskian. 



The same steps are followed to calculate x 2 for which we obtain 

i2(r) _ B(r) r «as<rv + KH /*' i™*- (46) 



In this case 



Y 1 (r)=r[K 1 (r)-^4§T I i(r)] 



h(R* ) 
Ki(R*) 
h(R*) 

g{r) = 2z 2 x„ (48) 



Y 2 (r)=r[K 1 (r)-^^h(r)} (47) 
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Having obtained the solutions x — xq + v 2 X2 and z = 1 + v 2 Z2 we can now consider the expression for the energy in 
this expansion: 



, 2 1 , I ^r/^O^ i ~2l i 4r (& Z 1\2 , 2,dx dX2 - 2 

£:(fi 2 )=7r?7 / dr <j — [( — ) +x ]+^[r( — ) +-( 



/ ;; I T dv 

.21 i / oC,.6\ 



dr rfr 



+ a; z 2 + X0X2) 



+ rz 2 2 ]+0(v 6 )} 
Using the equations of motion, the energy takes the form 



E(n) 

ITT] 2 



2 dxo 



(49) 



(50) 



The above expression gives us the energy of the vortex configuration up to the order v*. However an obvious 
drawback of this expression is the presence of the point Rq , which should be located numerically, and then preventing 
any analytical predictability power of the equation. Nevertheless we can convince ourselves that the point Rq can 
be shrunk to zero. The reason is that at leading order we can approximate the whole solution as a superposition of 
a Bogolmo'nyi vortex and a vortex concentrated at the boundary. Because both kind of solutions are exponentially 
small in complementary regions, any contribution from the non-linearity of the equations is exponentially suppressed. 
In fact, a simple plot of the solutions ( ^l| , ^4p^ ) shows that the solution with Rq — only differs in about 10 -3 with 
the one with non-zero Rq. Therefore, we will take Rq = 0. 

In this case the energy can be found by a simple numerical integration. Notice that the integrals only depends on 
R* and not on $ or n. The result can be compared with the obtained in reference 
a comparison between the energy values obtained from our approximate equation "(|50|), 

expression of reference and Q and the exact ones, corresponding to vortex solutions with R = 10 and n = 1. We see 
that even for big values of u, our approximate equation gives a result in excellent agreement with the exact ones. 



ice [El. In the following table we show 
tion (|q), the ones obtained from the 
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Energy 
from (50) 


Energy 
of Ref. j§ 


Exact 
energy 


1.0 


1.0 


1.0 


1.0 


3.0 


1.42 


1.4 


1.42 


5.0 


2.64 


2.6 


2.64 


8.0 


5.78 


5.9 


5.79 


12.0 


11.51 


13.1 


11.27 


18.0 


18.81 


29.9 


18.0 



Furthermore we can find a large R expansion of equation (|5(]). Using the asymptotic expansion of the Bessel 
functions, the first term in (50) has the asymptotic form v 2 (R* + 0.5 + 0(\/R*)). 



found by a numerical fit. The resulting approximate expression for the energy is 

E 

TTT] 

where 



, - n + v 2 R*a(R*) - v 4 R*b(R*) + 0{l/R\v 6 ) (51) 



a(R*) = 1 + l/(2R*) , b(R*) =0.139 + 0.111/E* (52) 

The most stable configuration corresponds to a vortex solution with vortex number n such that is an absolute 
minimum of the energy. Equation ( |5i"| ) can be approximately minimized with respect to n, with b(R*) / a(R*) as the 
expansion parameter. We find that the vortex number that minimizes the energy is given by: 
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R* 


( a 3 -b \ 


1" 


$0 


~ 2a 


Va 3 -3/26j 


+ 2 



(53) 



where [x] means the integer part of 




FIG. 3. The magnetization —M/(qr/ 2 ) as a function of the external flux $/<&o for R* — 10. 



The magnetization of the system shown in hg.3 is obtained from the Gibbs free energy, G = E — 2mf / Qq) 2 / R* 2 , 



as 



M 



dG tttj 2 ( 2a(R*) ( $ 



d<s> 



$0 \ R* 



4b(R*) ( $ 



$ 1 



(54) 



where n is given in equation (]53 

In this work we have analyzed the existence of vortex solutions of the Ginzburg Landau theory at the self dual 
point k 2 = 1/2 for a two dimensional disk of finite radius. Our original aim was to pursue further the interesting 
proposal made by Akkermans and Mallick Q of exploiting the properties of the BogomoPnyi equations for the study 
of the vortices at the self-dual point. Unfortunately our numerical study revealed that some of the assumptions made 
there are not entirely correct. 

We have shown that the minimal energy configurations do not satisfy the Bogomol'nyi equations in the inner disk. 
Nevertheless they provide a very good approximation to the actual solutions. Concerning the behaviour of the fields in 
the outer ring we have provided a simple analytical approximation scheme which do conform the numerical simulation 
and allows to obtain a selection rule for the number of vortices as a function of the external flux. 
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